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Abstract 


The  Fermi-Thomas  method  is  used  in  this  paper 
to  calculate  the  cohesive  energy  of  metals.  The 
preliminary  results  obtained  Indicate  that  there  is  a 
minimum  value  involved  in  the  relation  between  the 
energy  and  the  crystal  lattice  constant  of  a  metal. 

For  heavy  metals,  the  atomic  spacings  corresponding  to 
this  minimum  value  are  shown  to  be  close  to  the 
experimental  values;  however,  for  univalent  alkali  metals, 
the  theoretical' values  are  too  small,  and  the  oohesive 
energies  so  obtained  are  muoh  larger  than  the  experimental 
values.  The  author  tends  to  think  that  such  discrepancies 
are  due  to  errors  typioal  in  an  application  of  the  Fermi- 
Thomas  method  to  atoms,  and  that  they  are  not  particularly 
pertinent  to  the  application  to  metals. 


I.  Introduction 

Fermi  statistics  was  first  applied  to  atomic 
problems  by  Thomas^17  in  1926.  It  is  assumed  in  this 
theory  that  for  all  values  of  momentum  smaller  than  a 
fixed  p,  the  distribution  of  electrons  is  uniform  in  phase 
space,  l.e.,  in  each  element  in  phase  space  of  volume 
h 3  there  are  two  electrons;  further,  it  is  assumed  that 
there  exists  a  potential  energy  V  in  the  atom,  which 
depends  only  on  the  distance  r  from  the  nucleus,  and 
thus  by  means  of  the  Poisson  equation  one  oan  construct 
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the  proper  equation  for  the  potential  energy  in  the  atom. 
Practical  calculations  of  the  cohesive  energies  of  atoms 
were  curled  out  by  Milne  and  by  Baker ^->>  ;  their 
results,  are  greater  than  that  obtained  by  Hartree  s 
method  by  20  percent,  Corrections.for  boundary  effects 
and  ekcnahge  energy  made  by  Scott'  '  bring  the  calculat¬ 
ions  £lato  good  agreement  with  those  using  Hartree  s 
method  and  with  the  experimental  data;  the  error  involved 
is  less  than  3  percent.  The  success  of  the  application 
to  atoms  has  made  possible  the  application  to  metals. 

Since  in  the  application  of  the  method  of  Wigner  and 
Seitz (5 J  to  calculate  the  cohesive  energies  of  metals, 
it  is  possible  to  get  agreement  with  experimental  data 
only  for  alkali  metals,  and  the  application  to  multivalent 
metals  is  very  difficult,  one  would  think  that  the 
statistical  method  could  be  used  to  overcome  such  difficult¬ 
ies.  The  work  by  Slater  and  Krutter'  '  in  1935  showed 
little  prospect  of  obtaining  the  cohesion  of  metals; 
no  minimum  value  occured  in  their  curve  of  the  total 
energy  of  the  metal  versus  the  atomic  spacing;  the  results 
were  not  significantly  improved  even  after  correction  for 
exchange  energy.  Such  results,  according  to  them,  were 
due  to  the  omission  of  the  correlation  energy.  A  more 
systematic  study  begun  in  1936  by  Gombas '  shows  that 
the  problem  in  the  calculations  of  Slater  et  al,  lies 
not  in  the  omission  of  the  correlation  energy  but  in  the 
unjustified  treatment  in  their  statistics  of  attributing 
the  same  properties  to  the  valence  electron  and  the 
electrons  in  the  ions.  According  to  Gombas,  one  could  get 
results  in  agreement  with  experimental  data  and  the 
cohesion  of  the  metallic  lattice  points  if  one  regarded 
the  valence  electrons  as  free  electrons  distributed 
uniformly  in  metals  and  having  different  statistics  from 
that  of  the  electrons  in  the  ions,  and  took  into  account 
the  cohesive  energy  and  correlation  energy  of  free 
electrons,  the  mutual  interactions  of  ions,  together  with 
the  exchange  energy  of  the  valence  electrons  and  the 
electron  shells  of  the  ions,  etc.  Only  alkali  and  alkali- 
earth  metals  are  considered  in  Gombas  calculations;  the 
valence  electrons  of  such  metals  are  practically  all  free 
electrons,  and  so  the  results  are  close  to  those  obtained 
by  Har tree's  method.  However,  for  other  metals  such  an 
application  would  have  fundamental  difficulties  due  to 
the  lack  of  clear  distinction  between  valence  electrons 
and  the  electrons  in  the  ions.  One  finds,  for  Instance, 
in  the  application  of  the  Hartree  wave  function  to  copper, 
that  for  the  10  e  electronic  charge  of  the  d  shell  about 
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0.4  e  is  distributed  outside  the  atomic*  sphere  having  a 
volume  equivalent  to  that  of  one  lattice  unit.  In  other 
words,  the  d  shells  expand  into  comparatively  broad 
energy  bands  in  crystals,  and  the  d  electrons  would  tend 
to  play  a  role  almost, as  important  as  that  of 
the  valence  electrons'0'  as  far  as  cohesive  energy  is 
concerned.  There  is  no  sufficient  reason  to  make  a  sharp 
distinction  between  valence  electrons  and  inner-shell 
electrons.  Further,  results  relevant  to  practical  condit¬ 
ions  can  be  obtained  with  statistical  methods  only  when 
the  system  under  consideration  has  a  large  number  of 
particles;  and  it  is  hardly  meaningful  to  proceed 
statistically  with  only  a  handful  of  particles.  After  the 
separation,  made  by  Gombas,  of  the  valence  electrons  and 
inner-shell  electrons  into  different  categories,  it  is 
hard  to  understand  why  the  valence  electrons  as  individual 
particle  systems  could  obey  statistical  laws.  The  results 
of  Gombas  calculation,  though  achieving  some  quantitative 
agreement  with  experiments,  could  likely  be  a  sort  of 
accidental  coincidence,  void  of  firm  foundation. 

In  this  paper,  we  hope  to  find  a  general  method 
of  solving  the  problems  of  metals  of  many  valence 
electrons,  or  rare-earth  metals.  We  believe  that  the 
limitation  in  the  work  of  Slater  and  Krutter  is  due  neither 
to  the  omission  of  correlation  effects  nor  to  the  combined 
statistics  of  valence  electrons  and  inner-shell  electrons, 
but  that  it  arises  from  the  improper  boundary  conditions 
they  employed.  In  the  following,  we  will  try  to  set  up 
proper  boundary  conditions  to  solve  the  Ferml-Thomas 
equation,  and  obtain  better  results. 

II.  The  Ferml-Thomas  Equation  and  Its  Boundary 
Conditions 


In  the  application  of  the  sphere  approximation  method 
of  Wigner  and  Seitz,  we  set  up  a  sphere  of  radius  R  around 
every  atom,  determined  by  the  condition  that  the  volume  of 
the  sphere  be  equal  to  the  atomic  volume  and  that  every 
sphere  be  electrically  neutral,  with  the  mutual 
Interaction  between  spheres  neglected  in  the  calculations. 

We  need  only  to  calculate  the  energy  of  each  sphere.  From 
this  energy  we  subtract  the  energy  of  an  isolated  atom  and 

obtain  the  cohesive  energy. 

Suppose  in  every  atomic  sphere  we  can  define  a 
potential  function  V(rj  depending  on  the  distance  r  from 
the  nucleus  alone;  then  -eV(r)  is  just  the  sum  of  the 
potential  energies  of  an  electron  due  to  its  Interaction 


with  the  nucleus  and  with  other  electrons.  The  kinetic 
energy  of  the  electron  is  p2/2m  and  its  total  energy 
~@E  is  given  e imply  by 


2m*~  '  eT 


«eE 


In  the  classical  case  the  momentum  p  must  be  a  real 
number,  i.e. ,  -e(E-V )  must  be  a  positive  number,  in  order 
that  the  notion  be  well  defined.  The  total  energy  Is 
~eE0  corresponding  to  the  maximum  momentum  P  of  the  electron. 
Thus  -©Sq  is  Just  the  energy  of  an  electron  on  the  Fermi 
surface  in  metals.  From  the  above  formula, 

P  =  (E0  -  V)  J  *  . 


According  to  Fermi  statistics,  the  space  charge  density 
is 


p»  = 


:</i3 


("^E“2We<JP°"  V)j§* 


Substituting  in  the  Poisson  equation, 

V*F  .  -4npe. 

with  the  assumption  of  spherical  symmetry  of  the  electron 
distribution,  we  obtain 


XA.(r^)  - 

l2** 

A  [2 me(V  -  E0)V. 

(1) 

r*  dr  V  dr  / 

3 

ft* 

With 

1 

M  =  ( 

‘J>£\ 

*  0,4686  % 

» - r— '  A* 

<2) 

\ 

.1282/ 

* 

i 

11 

H 

1*0- 

<T  = 

if  =  2.2590  z*  ry. 

(3) 

*  a 

It  is  seen  from  the  above  formulas  that 


(V  «  Eo)r,  as  r  •#*  0,  r(V-Eo)  &e 


and  thus, 


C4) 


0  (0)  =  1. 

With  the  new  variable,  we  can  write  the  Poisson  equation 

&s  »  r 

'  "  X 


4>‘ 


<5) 
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According  to  Slater  and  Krutter,  the  condition 
that  the  atomic  sphere  he  neutral  simply  implies  that  at 
the  boundary  the  electric  field  vanishes#  Thus t 


dr  l»«ff 

where  R  is  the  radius  of  the  sphere,  That  la* 


=  ±*,  X  -  <«i> 

A*  P 

Here  ix  is  the  value  of  0  at  the  boundary  of  the  sphere. 
With  the  boundary  condition  (4),  one  can  then  proceed  to 
solve  Eq.  (5)  numerically,  (liven  an  initial  slope  » 

there  corresponds  a  graph  of  a  solution  p  versus  x. 

Fro®  condition  (6)  it  is  clear  that  the  tangent  to  the 
curve  mad©  from  the  origin  will  meet  the  curve  at  the  point 
whose  abscissa  corresponds  to  the  radius  of  the  atomic 
sphoi?©*  For  &  giv©3&  initial  slop©*  tSx©  solution  for  an 
isolated  atom  is  given  by  0  whose  graph  against  x  is 
asymptotically  tangent  to  the  x-axis  at  infinity.  This 
numerical  solution  provides  us  a  sequence  of  values  o* 

X,  0’o.  0*x  and  0x.  From  these  values  one  can  then 
calculate  the  potential  energy,  kinetic  energy  and  total 

energy  of  the  crystal.  .  . 

The  total  potential  energy  of  an  electron  in  the 

atom  is 

1  / '  P'Vdv  +  ±:eVn, 

2  Jo  2 


Pe  being  the  electron  density ,  and  Vo  being  the  electric 
potential  due  to  the  electron  at  the  position  of  the 
nucleus. 


pt 


i 

7  jj£ 

jpl' 


Va  = 


r  -  » 


I  i 

After  substitution,  with  help  of  $  * 

one  can  perform  an  integration  by  parts  similar  to  tha« 
of  Milne,  and  obtain 


Potential  , 

energy  „  j_  [* t  r  j£Vr£  +  Ea) ^  +  1 «  fx < - r >-£d*  - 

2  Jo  '  A*s  Jr'  X  f  2  j,; 
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_  -Iter  [*  **  *id*-Tz'rf'-iidl~ 

-  -}  [{ (**'-*<*'>’  + 1  **)j  - 

From  the  boundary  condition  0  (0)  a  1,  on®  can  obtain 

Potential  3  r  __  2.  xi  41  zer* .  < ' 7  > 

energy  3  VAJ  2  2 

The  maximum  kinetic  energy  of at3d 
e(V  ~  Eft).  The  mean  kinetic  energy  ^  » 

so  according  to  Fermi  statistics  the  kinetic  energy  i 

given  by  # 

-  [X  ~(V  -  E0)pedv  =  f  zer  ^  i  Mx  - 
Kinetic  energy  */o  5  J  *'  *a 

=  ~  zer  £-  <&  +  «].  * H 

Total  energy  »  potential  energy  +  kinetic  energy 

4]-izeF0~|zer«#.:v.  <») 

Mow  the  remaining  problem  is  to4.Je^®r®^?r5°g0  jS%U°ndlt 
ion  used  by  Slater  and  Krutter  to  determine  Eo  is  as 

SfXi.nT  Interatomic 

©pacing,  or,  the  electric  potential  in  the  J1®®0tljjlnlty 
of  a  nucleus  in  a  metal  should  be  equivalent  to  tha^ 
corresponding  to  a  nucleus  in  an  Isolated  atom.  Thus 

obtain 

E0  *  rW«°)  -  ^o<X)l. 

Hera  £«*»  is  the  initial  slope  corresponding  to  the 

solution  for  an  isolated  atom;  and 


£zdv  —  — r4te<°°>* 
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It  follows  that 


Total  energy  =  ~  z<t 


\4>' 


15 


*  4 

i)' 


(10) 


Starting  with  this  formula  '  one  falls  to  obtain 
a  minimum  value  for  the  curve  of  the  total  energy  against 
X;  In  other  words,  fro®  the  above  on©  is  not  able  to 
illustrate  the  stability  of  the  cohesion  of  metals . 

In  our  opinion  the  above  boundary  condition  is  not 
adequate.  The  electric  potential  of  an  electron  is  the 
near  neighborhood  of  the  nucleus  is  the  work  done  In 
taking  an  electron  from  infinity  to  that  neighborhood. 

This  should  obviously  be  different  for  an  isolated  atom  and 
for  an  atom  in  a  metal.  To  bring  an  electron  from  infinity 


to  the  inside  of  a  metal  on©  has  to  overcome  the  boundary 
effect;  in  other  words,  the  actual  work  done  is  different 
from  that  for  an  Isolated  atom  by  the  work  function  of 
the  metal.  Slater  ©t  al.  overlooked  this  point,  and  so 
could  not  obtain  a  stable  metallic  cohesion. 

After  some  calculations,  we  have  reached  the 


conclusion  that  their  assumption  together  with  the  assumpt¬ 
ion  of  zero  potential  at  the  boundary  of  the  sphere  is 
exactly  equivalent  to  the  following  assumption: 


Vo  —  fx**v~  ~T(<t>x  -  &> 

r 

the  latter  assumption  leads  to  results  identical  to 
Eq.  (10).  However,  this  latter  assumption  was  taken  by 
slater  et  al.  as  being  of  little  physical  meaning. 

Fermi -Thomas  statistics  takes  into  consideration 
all  the  electrons  in  its  calculation  of  electric  charge 
density  at  a  point  in  space.  When  we  calculate  the 
electron  potential  energy  from  this  charge  density,  we 
should  notice  that  the  potential  energy  is  the  energy 
of  the  electron  due  to  the  potential  field  of  the  nucleus 
and  other  electrons.  Therefore,  rigorously  speaking,  to 
find  the  potential  energy  of  a  given  electron,  one  should 
proceed  with  a  statistics  that  considers  the  potential 
energy  of  all  electrons  except  for  the  given  one. 

However,  to  find  the  electric  charge  distribution  one 


lT~ (10)  seems  to  involve  some  error.  Since  in  Slater's 
paper  the  actual  calculations  are  not  fully  presented, 
we  cannot  say  whether  Slater's  calculation  is  accurate 
or  not. 
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has  to  consider  the  effect  of  this  particular  electron  on 
all  other  electrons.  In  this  manner,  the  potential  energy 
of  the  given  electron  can  be  found  from  the  electric  charge 
density  distribution.  Since  our  statistical  object  is  a 
large  system  of  2  electrons,  the  density  distribution  of 
2  -  1  electrons  can  be  substituted  by  that  of  2  electrons} 
the  resulting  fractional  error  in  the  substitution  is  not 
larger  than-*-.  But  in  considering  the  electric  potential 
of  a  given  *  electron  at  the  boundary  of  the  atomic  e 
sphere  due  to  all  other  electrons ,  one  should  take 
and  not  zero.  Since  the  potential  seen  by  the 
given  electron  is  due  to  2  -  1  negative  and  2  positive, 
charges,  the  electric  potential  at  the  boundary  cannot  be 
regarded  as  vanishing  as  in  the  case  of  a  potential  of  an 
averagely  neutral  sphere. 

The  boundary  condition  is  thus: 


_  e  „  JL  =  r 
Tssti  — rr  v  1  ~  V 

R  g  A 


So 


After  substitution  in  Eq.  (9)  and  some  manipulation  we 
obtain 

Total  energy  -  set  [*S  +  —  X*  £  ~  j yj-  ( 5 


III.  Results  of  the  Calculation  _of  thg. 

Cohesive  Energy 

The  cohesive  energy  of  a  metal  is  the  difference 
between  the  total  energy  of  an  atom  in  the  metal  and  the 
total@energy  of  the  isolated  atom.  The ' ordgr  of  magnitude 
of  the  cohesive  energy  is  about  10  to  10  J^at  of 
the  total  energy  of  an  isolated  atom.  Therefore,  in  older 
to  calculate  accurately  the  cohesive  energy  one  J*®,  J*®.  t 
solve  Eq.  (5)  accurately  so  as  to  get  about  eight  significa 
figures.  So  far  there  is  no  such  accurate  numerical  ^ 
solution.  In  the  following  we  will  employ  two  methods  oa 
calculation  and  compare  the  results.  These  two  methods  are 
Stonily  to.  same .  We  shall  first  eliminate  the  energy  of 
the  isolated  atom  by  an  approximate  expansion*  ana 
then  correct  our  calculation  with  help  of  the  numeriea 
solution  of  Eq.  (5). 
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1.  Approximate  expansion 
Set  $  **  fa  4>» 

Here  0*  is  the  function  for  the  isolated  atom, 
is  an  additional  function.  They  are  of  the  same  order  of 
magnitude  for  large  r.  However,  we  take  9 
general  and  obtain  the  approximate  solution  oy  substitut¬ 
ing  the  above  relation  in  Eq.  (5)t 

X 1 

In  the  isolated  atom,  0a  satisfies  the  equation: 


*A 


A 


*• 


(12) 


and  hence, 


0" 


(*)** 


(13) 


One  can  make. use  of  the  solution  to  Eq.  (12)  by 
Rush  and  Caldwell (^^  to  solve  Eq.  (13)  numerically. 
Equation  (13)  is  a  linear  differential  equation,  and  so 
Its  solution  may  be  different  by  an  arbitrary. oon|^tt® 
which  can  be  determined  by  boundary  conditions.  From  th 
boundary  condition. 


0x 


i*.  <t>'A  + 


4>Ar  +  B<ft, 


we  obtain 


B 


0AX  — 

X</4  -  0.V 


(14) 


B  is  a  function  of  the  radius  X  of  th® 

On  substituting  the  above  expressions  in  Eq.  Ul J 
eliminating  the  function  for  Isolated atom,  we  obtain  the 
cohesive  energy  AE  in  terms  of  the  atomic  radius 

X: 

A  £  -  A  zer  [fitf  +  £x*(4>a  +  B4>xy  -  -f 

7  L  15  *  6  ZAJ 
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There  Is  a  minimum  value  in  the  graph  of  .*f*inst  X* 
The  radii  of  the  atomic  spheres  for  \ 
corresponding  to  this  minimum  value,  are  shown  in  -aui . 

Table  1  (R  in  units  of  &  ) 


Na  j  A1 
ii  !  u 


Cr 

Co 

Cu 

"  Mo 

Ag 

T& 

Pt 

24 

2? 

29 

42  ! 

47 

73 

78 

1.41 

1.39 

,.41 

,.54 } 

1.S9 

i  1*1 

i 

1.S4 

- - - 

j  1.39 

V35 

1.38 

1.3*  i 

1.32 

j  1.26 

1.23 

Element 

Atomic  number  i  ,  n 
Experiment  |  2 0 
K  Theory  r*  . 


The  values  of  EE  corresponding  ^oi-heminlmum 
-»  __  oil  fno  toY  1  OX*  2  OX%d©X'S  of  ffiEgllltlUa  • 

Sis  is  due  to  ?he  approximation  in  the  calculation  for 
theS Isolated  atom?  F?om  the  nature  of  the  approximation 
one  can  roughly  estimate  such  errors. 

2.  Numerical  evaluation 

Se"^oiran“frhM,\  Maroh(H]  and  by  Latter^, 

Jia Among  these  the  calculation  by  Latter  is  most 

fo?  devaluation  of  the  cohesive  energy.  Here  we  have 
St*!  tables  to  find  the  relation  between  0o, 
f*  Al  l  Jhile  we  estimate  two  more  figures.  The  numerical 
SsSns  are  as  sSown  in  Table  2.  The  initial  slope  is 
about  -  1.588073  for  the  Isolated  atom  solution. 
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Table  2 


X 

Qx 

X 

-e; 

4.905 

1.58376ft 

0.20099 

i  *.733 

1.587854 

0.07!  m 

$.229 

1.584923 

0.18058 

9.565 

1.587932 

0059868 

ft  0 1 

1.585302 

0.17095 

9.858 

1.587970 

0.056426 

5.506 

t. 585504 

0. 1651 1 

10.804 

1.588023 

0.046934 

S.850 

1.586275 

0.14899 

11.300 

1.588030 

0.042832 

6177 

1.586675 

0.13547 

11.963 

1.588033 

0.038086 

6.20* 

1.586687 

0.13433 

12.369 

1.583036 

0,035529 

7.014 

1.587314 

0.10539 

13.771 

1.588040 

0.028322 

7.385 

1.587435 

0.098098 

15.471 

1.588045 

0.02202! 

7.790 

1.5*7560 

0.088823 

15.870 

1.588050 

0.020825 

8.015 

1.587643 

0.084193 

16.00 

1.588060 

0.020454 

8.588 

1.587798 

0.073821 

If  we  substitute  in  Eq.  (11)  the  numbers  given  in 
Table  2,  we  can  obtain  a  graph  of  AE  against  X;  this 
graph  has  a  minimum  value  which  corresponds  to  atomic 
radii  as  shown  in  Table  3. 


Table  3  (R  in  units 


of 


i) 


Element 
Atomic  number 
Experiment 
R  Theory 


|  Na 

j  A1  ! 

K 

V 

Cr  | 

Co 

Cu 

'  Mo 

|  Ag 

i  •» 

|  ,9 

23 

24  | 

1 

27 

I  29  | 

42 

47 

j  2.08 

j  1.62 

I _ _ _ _ 

j  2.57 

j  1.47 

1.41  ! 

S  1.39 

i  i.4i 

1.54 

1.59 

1  1.47 

1 

{  1.49 

|  1.41 

j  1.37 

1.38 

1.37 

M0 

1.52 

1,57 

The  A/;  values  obtained  are  equivalent  to  those 
obtained  by  the  previous  approximation  method. 


IV.  Discussion 


Prom  the  above  tables  we  see  that  the  atomic  radii 
derived  from  the  theory  are  basically  in  agreement 
with  the  experiments.  For  multivalent  and  B-group  metals, 
e.g. ,  V,  Or,  Co,  Cu,  Ag,  etc.,  the  agreement  is  very  good, 
but  for  A-group  univalent  alkali  taetals,  Na,  etc**  the 
deviation  becomes  larger.  Since  the  cohesive  energy  of 
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alkali  metals  Is  due  chiefly  to  the  change  In  electron 
cloud  distribution  of  one  single  valence  electron,  the 
reason  for  the  deviation  is  easily  .seen  to  be  the  large 
error  caused  by  the  application  of  Fermi  statistics  to 
the  single-electron  system.  Therefore  one  should  not 
expect  consistent  results  for  univalent  alkali  metals. 

As  to  the  great  deviation  between  the  theoretical 
and  experimental  vahes  of  cohesive  energy.  We  would 
attribute  It  to  the  following  reasons  the  cohesive 
energy  of  metals  is  of  the  order  of  10"4  to  10"6  that  of 
the  total  electronic  binding  energy  of  the  atom;  a 
1  percent  or  even  0.1  percent  error  in  the  total  energy 
of  the  atom  is  enough  to  cause  a  completely  erroneous 
result  in  the  cohesive  energy.  Since  the  '  Fermi-Thomas 
method  is  accurate  to  the  order  of  ,  it  would  be  hard 
to  do  any  accurate  numerical  calculation  of  the  cohesive 
energy.  If  one  can  find  a  mathematical  technique  for 
cancelling  the  energy  of  the  atom  from  the  final  result 
so  that  it  will  not  appear  in  the  final  result  for  the 
cohesive  energy,  then  one  can  eliminate  a  great  portion 
of  the  error. 

In  the  Fermi-Thomas  statistical  method,  we  use  the 
classical  concept  that  the  electron  momentum  p  has  to  be 
real  so  that  the  electron  kinetic  energy  will  be  positive. 
This,  however,  is  not  consistent  with  the  quantum- 
mechanical  concept.  Quantum  mechanics  Indicates  that  in 
the  outer  part  of  the  atom,  the  kinetleenergy  of  the 
electron  can  be  negative.  When  the  Fermi-Thomas  method 
is  applied  to  the  calculation  of  the  total  energy  of  the 
atom,  no  significant  error  could  occur  from  omitting  this 
negative  part  of  the  kinetic  energy;  however,  the 
numerical  value  of  the  cohesive  energy  could  be  affected 
immensely.  Therefore,  in  an  accurate  calculation, 
corrections  should  be  made  according  to  an  estimate  of 
magnitude  of  this  negative  kinetic  energy.  In  actual 
problems,  for  example,  the  calculation  of  energy  of  an 
alloy,  one  can  make  use  of  the  experimental  data  for  pure 
metals  to  decide  on  this  part  of  the  correction. 

This  paper  presents  a  preliminary  calculation  of 
the  cohesive  energy  of  pure  metals  without  taking  into 
consideration  the  exchange  energy  and  the  correlation 
energy.  Nevertheless,  the  results  obtained  prove  to  be 
an  improvement  over  the  work  by  Slater  and  Kr utter,  since 
the  negative  minimum  value  for  the  cohesive  energy  of 
metals  and  the  radii  of  the  atomic  spheres  obtained  are 
basically  in  agreement  with  the  experimental  data.  It 
Is  hoped  that,  after  some  modification  of  the  above 
method,  it  would  be  possible  to  calculate  the  cohesive 
energy  of  heavy  metals,  e.g. ,  the  rare-earth  metal  alloys 
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and  the  hlack  metal  alloys.  It  is  therefore  of  considerable 
significance  for  the  study  of  alloys. 
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